Vacancy-induced spin textures and their interactions in a classical spin liquid 
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Motivated by experiments on the archetypal frustrated magnet SrCrg p Gai2-9 P Oi9 (SCGO), we 
study the classical Heisenberg model on the pyrochlore slab (Kagome bilayer) lattice with site- 
dilution x = 1 — p. This allows us to address generic aspects of the physics of non- magnetic vacancies 
in a classical spin liquid. We explicitly demonstrate that the pure (x = 0) system remains a spin- 
liquid down to the lowest temperatures, with an unusual non-monotonic temperature dependence 
of the susceptibility, which even turns diamagnetic for the apical spins between the two kagome 
layers. For x > but small, the low temperature magnetic response of the system is most naturally 
described in terms of the properties of spatially extended spin textures that cloak an "orphan" 
S — 3/2 Cr 3+ spin in direct proximity to a pair of missing sites belonging to the same triangular 
simplex. In the T — > limit, these orphan-texture complexes each carry a net magnetization that is 
exactly half the magnetic moment of an individual spin of the undiluted system. Furthermore, we 
demonstrate that they interact via an entropic temperature dependent pair-wise exchange interaction 
J e f f (T, r) ~ TJ{fy/T) that has a logarithmic form at short-distances and decays exponentially 
beyond a thermal correlation length f(T) ~ 1/VT. The sign of J e // depends on whether the two 
orphan spins belong to the same Kagome layer or not. We provide a detailed analytical account 
of these properties using an effective field theory approach specifically tailored for the problem at 
hand. These results are in quantitative agreement with large-scale Monte Carlo numerics. 

PACS numbers: 75.10.Jm 05.30.Jp 71.27. +a 



I. INTRODUCTION 

Magnetic ions in Mott insulators often interact with 
short-ranged antiferromagnetic exchange couplings J. 
When ions with spin S occupy a bipartite lattice, the 
system develops antiferromagnetic order with the spins 
forming a collinear Neel state below a transition temper- 
ature Tjv of order the Curie- Weiss temperature @cw ~ 
JS 2 . However, if magnetic lattice defined by the nearest 
neighbour connectivity matrix of the magnetic ions has 
triangular loops (more generally, loops with an odd num- 
ber of sites) , the antiferromagnetism is frustrated, and the 
system does not develop long-range Neel order. 

In such cases, the classical exchange energy typically 
possesses a very large number of minima — indeed, the 
ensemble defined by the classical minimum- energy con- 
figurations often has finite entropy in the thermodynamic 
limit. This leads to "spin-liquid" behaviour over a broad 
temperature range T/ <Tc @cw in which the physical 
properties reflect averages over all possible minimum en- 
ergy configurations, and are therefore largely determined 
by the geometry of the magnetic lattice. Below the freez- 
ing temperature Tf, the system eventually orders, but 
the ordering patterns are often quite complex and deter- 
mined by the interplay between subleading terms in the 
interaction Hamiltonian and the effects of quantum and 
thermal fluctuations.— 

In the spin-liquid regime, the geometric frustration ef- 
fectively "quenches" the leading exchange interactions — 
in some ways, this is similar to fractional quantum hall 
systems in which the formation of Landau levels quenches 
the leading kinetic energy term in the Hamiltonian. As in 



the fractional quantum hall case, this can lead to unusual, 
emergent degrees of freedom dominating the physical re- 
sponse of the system^ One well-studied example of this 
is the low temperature physics of spin-ice compounds in 
the classical regime, which admits a natural description 
in terms of emergent magnetic monopole degrees of free- 
dom in a classical easy-axis magnet^ Another quantum- 
mechanical realization of such quasiparticle fractionaliza- 
tion is the Coulomb liquid phased— e.g. of bosonic mat- 
ter on the three-dimensional pyrochlore lattice. 

Here, we focus on the physics of the frustrated an- 
tiferromagnet SrCr 9j ,Gai2-9pOig (SCGO)jSr— in which 
the Cr 3+ S = 3/2 ions interact with nearest neighbour 
Heisenberg exchange interations J « &0K to form a 
corner-sharing network of antifcrromagnctically coupled 
spins that can be variously described as a Kagome bi- 
layer or a pyrochlore slab lattice whose sites are diluted 
with a density x = 1 — p of vacancies; these vacancies 
reflect the presence of non-magnetic Ga ions on the Cr 
sites in SCGO for p < 1. Using the classical Heisen- 
berg antiferromagnet on this SCGO lattice to model the 
low temperature physics, we demonstrate that the ex- 
perimentally hitherto unachievable pure compound with 
p = 1 remains a spin-liquid down to the lowest temper- 
atures accessible in classical Monte-Carlo (MC) simula- 
tions; this is in keeping with analytical expectations.— 

We also demonstrate that the low temperature mag- 
netic response of the site-diluted system with dilution 
x is most naturally described in terms of the proper- 
ties of spatially extended spin textures that cloak "or- 
phan" 5 = 3/2 Cr 3+ spins in direct proximity to a pair 
of missing sites belonging to the same triangular sim- 
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plex. Our Monte- Carlo studies establish that a single 
orphan-texture complex in the low temperature limit has 
a net moment equal to exactly half the moment of a 
free S = 3/2 spin. Roughly speaking, this arises from 
the fact that the orphan spin "sees" an effective mag- 
netic field h/2 upon application of a uniform external 
field h — half the external field being "screened" by the 
exchange-coupling of the orphan spin to the surrounding 
spin-liquid — and the surrounding spin texture "cancels 
off" exactly half of the orphan spin's magnetization by 
developing a net diamagnetic response of the right mag- 
nitude in the low temperature limit. 

This "fractional moment" is in agreement with predic- 
tions of an effective field theory approach that we de- 
velop here. Furthermore, as already noted in our earlier 
Letter—, the asymptotic low temperature behaviour pre- 
dicted by this effective theory is found to be surprisingly 
robust in MC simulations, persisting up to temperatures 
of order 0.1J5 2 for spin-S 1 magnets on the SCGO lat- 
tice with nearest neighbour coupling J. Here, we show 
that the underlying reason for this robustness has to do 
with the fact that this fractional moment of an individ- 
ual orphan-texture complex at T = is, in a well-defined 
sense, much more localized than the l/|r| far-field fall-off 
of the surrounding spin texture. 

At finite density, these textures (here and henceforth, 
we slur over the distinction between the orphan-texture 
complex and the texture unless absolutely essential) in- 
teract with each other by an emergent temperature depen- 
dent pair-wise exchange interaction J e j / (T, r) that de- 
cays exponentially beyond a thermal correlation length 
£(T) - 1/VT and whose sign depends on whether the 
two orphan spins belong to the same Kagome layer or 
not. Again, the absence of three-body and higher inter- 
actions, as well as the dependence of J e // on T, r and 
layer index seen in MC simulations are all successfully 
modeled using the effective field theory approach we de- 
velop here. 

Thus the physics of such vacancy-pairs dominates the 
impurity contribution to the susceptibility of the Kagome 
bilayer at small x in the low temperature limit. Although 
the contribution of these orphan spin textures is masked 
in macroscopic susceptibility measurements by another 
"extrinsic" contribution (see Sec. [TTJ) that has little to do 
with the physics of the frustrated Kagome bilayer, these 
spin textures play a crucial role in determining Knight 
shifts and lineshapes in Ga NMR, as has already been 
emphasized in our previous Letter 

This article is organized as follows. Sec. [TT| provides a 
whistle-stop tour of the most pertinent previous theoret- 
ical work. Our main technical tools, including the effec- 
tive field theory approach we develop here, are described 
in some detail in Sec. IIIII and Sec. IIV1 which may be 
skipped by a reader only interested in final results of rel- 
evance to SCGO. The following sections IVl IVII contain a 
description of such results, first for the pure system, then 
for a single texture, and finally for interactions between 
the textures. The article closes with considerations on 



the generality of our results — in particular, how to extend 
them to other lattice topologies and dimensionalities — 
along with an outlook on future work. 



II. BACKGROUND: CLASSICAL HEISENBERG 
MODEL ON THE SCGO LATTICE 

We consider classical length-5 1 spins interacting with 
nearest neighbor Heisenberg exchange couplings on the 
SCGO lattice (see Fig [J): 

H nn = J Sf ■ — h ■ Sft (f ) 

(rr) r 

Here J > is the antiferromagnetic Heisenberg exchange 
interaction on a link (rr ) connecting nearest neighbour 

sites f and r on the SCGO lattice, and h is the uniform 
external magnetic field. 

As depicted in Fig. [TJ the Cr 3+ ions in SCGO define 
a lattice consisting of a pyrochlore slab with one layer 
of up-pointing tetrahedra sharing a vertex each with a 
second layer of down-pointing tetrahedra. Thus, this 
Kagome bilayer or pyrochlore slab may be thought of as a 
corner sharing network of up-pointing tetrahedra, down- 
pointing tetrahedra, and triangles that do not belong to 
any tetrahedron — the latter are the would-be bases of 
the next layers of tetrahedra which would be needed to 
convert the SCGO lattice into a three dimensional py- 
rochlore structure. These next layers of tetrahedra are 
absent in the SCGO structure, being replaced in effect by 
a layer of magnetically "inert" Cr dimcrs (spin pairs) that 
serve to isolate each Kagome bilayer from the Kagome bi- 
layer in the next unit cell, and lead to an effectively two 
dimensional situation. The main effect of this layer of 
Cr dimers is to give rise to a population of isolated Cr 
spin S = 3/2 moments with density proportional to x 
when Ga impurities replace an 0(x) fraction of Cr spins 
in this isolated dimer layer. Although the resulting pop- 
ulation of free Cr spins gives a dominant contribution to 
the macroscopic susceptibility at small x and low tem- 
perature, this has nothing to do with the physics of the 
frustrated Kagome bilayer. Therefore, we leave the Cr 
dimer layer out of our discussion in the rest of this arti- 
cle. 

Keeping this in mind, we rewrite the classical Hamil- 
tonian of the system as 

ms}) = { BE & - ijf + { »E & - y 

GS rGKI A re A 

(2) 

where H and A denote the basic tetrahedral and trian- 
gular simplical units of the kagome bilayers in the SCGO 
lattice (Fig[T]). It is now clear that this classical exchange 
energy has a global minimum when each simplex of the 
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for each spin component a. At zero magnetic field, this 
implies that the total spin of each corner-sharing simplex 
is zero. 

The T = 0, h = physics is dominated by the en- 
semble of states which achieve this global minimization 
of the energy by setting the total spin of each simplex 
to zero, and this is expected^— to be the case even in 
the presence of site dilution: this ansatz of "simplex sat- 
isfaction" , whereby the total spin of each simplex is set 
to zero in the configurations that dominate the T = 
physics, has dramatic consequences. If any simplex has 
more than one spin remaining on it after dilution, the 
corresponding simplex continues to have zero total spin, 
while a "defective simplex" with only one surviving spin 
on it cannot have zero total spin. As a result, an infinites- 
imal magnetic field acting at T = in the z direction only 
couples to such defective simplices, yielding the follow- 
ing saturation magnetization for a system with a single 
defective simplex: S? ot = ±(£ K 5| + £a S' a ) = S/2. 

Thus, a single defective simplex leads to a magne- 
tization S/2 in response to an infinitesimal magnetic 
field at T = 0, hinting at the presence of fractional- 
ized spin degrees of freedom liberated by such vacancy 
configurations. Although this argument is explicitly re- 
stricted to T = 0, and although it is clear that spin- 
correlations at low temperature extend over large dis- 
tances in SCG0 16 i 17 , it is curious to note that precisely 
these defective simplices give rise to a low temperature 
Curie response cx 1/T within the "single- unit approxima- 
tion" of Moessner and Berlinsky that does not explicitly 
incorporate spin-correlations beyond a single simplex. 

Is this a coincidence, or are fractionalized spin S/2 de- 
grees of freedom "real" in the low temperature classical 
spin-liquid regime in SCGO and related systems? Here, 
we address this question using classical Monte-Carlo sim- 
ulations of the microscopic system as well as an effective 
field theory approach that correctly incorporates entropic 
effects on the same footing as the energetics of simplex 
satisfaction. We will see that this approach enables us to 
describe the physics of vacancies with a remarkable level 
of accuracy and analytical detail. 



III. METHODS: EFFECTIVE FIELD THEORY 
AND CLASSICAL MONTE-CARLO 

Here, we provide an overview of the effective field the- 
ory and classical Monte-Carlo methods we use to explore 
the physics of orphan spins in SCGO. The reader not 
interested in technical details may skip this section. 
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FIG. 1. The Cr atoms in SCGO form a lattice of Kagome bi- 
layers separated from each other by a layer of isolated dimers 
consisting of pairs of Cr atoms. Each Kagome bilayer is a cor- 
ner sharing arrangement of tetrahedra and triangles made up 
of two Kagome lattices coupled to each other through "apical" 
sites shared between up-pointing and down-pointing tetrahe- 
dra. Links between near-neighbour sites in each bilayer rep- 
resent a Heisenberg exchange coupling J = 80 K between 
neighbouring Cr 3+ spins, while links in the isolated dimer 
layer represent a Heisenberg exchange coupling J' — 216 K 
between the two Cr 3+ ions that constitute each pair. Two va- 
cancies in the (shaded) triangle leave behind an orphan spin. 



A. Monte Carlo simulations 

In order to test our theoretical predictions, we employ 
detailed classical Monte Carlo simulations to study Eq. Q] 
with 0(3) and 0(4) spins. To update a spin configura- 
tion efficiently, we use^S a combination of "microcanon- 
ical" moves that conserve energy and the usual "canon- 
ical" Monte-Carlo moves. In all cases, we use periodic 
boundary conditions, and study a sequence of sizes rang- 
ing from L = 16 to L = 50 (with 7L 2 sites) to control 
finite-size effects and reliably extrapolate to the large L 
limit. In runs that focus on autocorrelation properties, 
we are careful to "switch off" the microcanonical moves 
to guard against the possibility that they mask glassy 
slowing down. 



B. Field theory: large-iV and soft spins 

In the pure case, the basic idea of the effective field 
theory is to treat the fixed length constraint Si = S 2 in 
an approximate way by replacing the original fixed length 
spins by fields <j)p whose average length is controlled by 
a phenomenological stiffness parameter to ensure that 
(</>!) = 5 2 . Thus, the partition function Z c ff of the 
effective theory for the pure system is given by 

Z cS = J T>$cxp(~S cff ) 

s c « = \Y,Pr$l + ^n ({&}) 

r 

(4) 
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where V(j> = ]J p a d(fip , and the stiffness constants are 
fixed as pp = pi (pp = P2) for ali the sites on the kagome 
(apical) layers by requiring that (</> 2 )z off = S 2 . Solving 
these constraints in the thermodynamic limit, we obtain 
Pl w l^SSOS" 2 and p 2 « 1.5903S 2 at T = and h = 0. In 
the spin liquid regime (T <C JS 2 , /i <C JS), the changes 
to /?! and P2 induced by non-zero temperatures and fields 
are small and can be ignored to a good approximation^. 

Although our vector field 4>p is a three component ob- 
ject defined at each lattice site, this effective field the- 
ory treatment of the pure system is equivalent to the 
leading term in a large- N expansion for the O(N) gen- 
eralization of the original classical Heisenberg spin sys- 
tem. Our rationale for relying on this approach is that 
this leading term in the large-iV expansion is known 
to provide a very good approximation to the physics 
of the classical Heisenberg antiferromagnet on both the 
Kagome and pyrochlore lattices in their respective spin 
liquid regimes i 18 ' 19 Since the SCGO lattice is essentially 
a finite-thickness slab of the pyrochlore lattice, we expect 
that this large- TV treatment will be a good approximation 
in the spin liquid regime on the SCGO lattice as well. In 
what follows, we explicitly test this assumption against 
the results of classical Monte Carlo simulations before 
proceeding to study the effect of vacancies. 

The fact that SCGO is a spin-liquid in this approxi- 
mation can be seen by taking the T — > of the action in 
Eq |U Then, we get 

Z eff = f Vfoexp \-\Y,P?%) II^ A *)II 

J V f / Af. § ? 

(5) 

The local constraints on the ground state manifold 
(<?W = EfeA = °> and fc ? = SfeKl = °) can 
be elegantly described by constructing vector fields B a 
that live on the bonds of the dual lattice of the SCGO 
lattice; the sites of this dual lattice are located at the cen- 
ters of the simplices of the SCGO lattice, and the links 
of this dual lattice pass through the SCGO sites. Thus 
we write 

B a (r) = e?<p a (f f ), a = x,y,z (6) 

where the unit vectors e are always directed from the 
A sublattice site of the (bipartite) dual lattice to the 
neighbouring B sublattice site. With this choice, these 
constraints translate to the requirement that each of the 
B a fields are divergence-free: 

V • B a = (7) 

and the action in Eq [5] becomes 

l^pAB?) 3 (8) 

T 

This implies spin correlators which are dipolar and decay 
as 1/r 2 in real space. The zero-temperature formulation 



and behaviour of such a Coulomb phase is well under- 
stood in related settingSi 16 i 17 ' 19 At non-zero but small 
temperatures, there is a finite thermal correlation length 
£ ~ 1/ y/T beyond which the spin correlations decay ex- 
ponentially In Coulomb language, 0a and <f>^ can be 
thought of as "vector-charges" , and these are strictly zero 
at T — but can be generated thermally, without an 
excitation gap, when T / 0. As we will see, the spin- 
"charge" correlator of the undiluted parent spin liquid 
determines the texture induced around an orphan spin 
and the "charge" - "charge" correlator determines the in- 
teraction between orphan spins within our effective the- 
ory. 



C. Field theory: Vacancy effects 

To handle the presence of vacancies, we employ La- 
grange multipliers A", i.e. the Fourier representation 

S(^) = ^J dXpexp(i\p^) (9) 

to enforce the constraint </>2 = for each component a 
of 4>p at all sites r = r v on which the Cr 3+ spins have 
been substituted by non-magnetic Ga. Beyond this, our 
central technical innovation lies in a careful treatment 
of those defective simplices, in which all spins but one 
have been substituted by non-magnetic Ga. It is abun- 
dantly clear that the fixed length constraint on the mi- 
croscopic spins is crucial for the single "orphan" spin left 
behind at the neighbouring site f on such defective sim- 
plices. In order correctly to capture the physics of these 
orphan spins, we retain the original microscopic fixed 
length spin-S* variable Sp o at each such orphan sites r a , 
i.e. we do not approximate these orphan spins by Gaus- 
sian fields <j)p o . 

Thinking in terms of the large-W limit of the O(N) 
spin model with vacancies provides an alternate perspec- 
tive on our approach: In the presence of vacancies, there 
is no translational invariance, and the stiffness parame- 
ters pp of the large- N theory acquire a dependence on the 
spatial position f. Since the most dramatic spatial de- 
pendence of pp is expected to arise at and perhaps around 
each orphan site r = r D , retaining these orphan spins as 
microscopic fixed- length vectors is a way of incorporating 
the leading effects of this spatial inhomogeneity. Having 
done this, we expect that the stiffness pp does not devi- 
ate significantly from its uniform system value at other 
sites. With this in mind, we ignore the r dependence of 
pp once the orphan spins are retained as spin-S" vectors. 
The validity of this assumptions is confirmed a posteriori 
against numerics. 

Our effective theory thus reduces to a field theory with 
a constrained Gaussian field <j)p coupled to fixed-length 
spin-S" vectors at the orphan sites. We are thus led to a 
sort of "classical Kondo lattice model" in which one has 
spin-S variables coupled to a bosonic "bath" represented 
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by 0. It is interesting to note that the effective coupling 
between the orphan spins, and the effective "g-factor" 
with which the orphan spins couple to the external mag- 
netic field, are both determined by the coupling of orphan 
spins to the bosonic field 0, which may be thought of as 
representing the degrees of freedom of the bulk spin liq- 
uid. 

Calculations within this framework are most conve- 
niently performed by noting that the fixed-length spin-S* 
character of the orphan spins can also be implemented 
by introducing additional Lagrange multipliers p~ to en- 
force the constraint (j)" = Sri- , where n? is a unit vector 
at each orphan spin site r Q (so that Sn? o is the length-.? 
orphan spin vector). Note that these Lagrange multipli- 
ers A and p play a quite different role in our theory from 
the stiffness parameters p: The p capture the effects of 
the fixed length constraint on average and are not inte- 
grated over, while the Lagrange multipliers A and p are 
integrated over in order to set cj) to at vacancy sites and 
Sn at the orphan sites. 

To avoid cluttering notation with factors of J and S, 
we now measure temperature in units of JS 2 (i.e. replace 
T/JS 2 by T) and magnetic field h a in units of JS (i.e. 
replace h a /JS by h a ) and use (4> 2 ) = 1. In other words, 
we write 



Vn J VX J Vp J V$exp(-S' cS ) 



S, 



eff 



S, 



eff 



where VX = H^dX^, Vp = Y\r^ a d ^h and Vri = 
Y\p a dnp S(n 2 —1). Since all these Lagrange multipliers 

\ft v and pf o couple linearly to the Gaussian field <f>, and 
since <S e fj({</y}, h, T) is the action for the translationally 
invariant pure system, one can do the Gaussian integrals 
over the 4>r fields exactly by performing the correspond- 
ing integrals over the Fourier modes of the <j) field. 



(10) 



Upon doing these integrals, one obtains an action writ- 
ten as a quadratic form in terms of the uniform external 
magnetic field h and the Lagrange multiplier fields A" 
and pp at vacancy sites r v and orphan sites r respec- 
tively. This quadratic form is of course diagonal in the 
spin index a, and h only appears in the action for the 
z component (a = z) if the field acts along the z axis. 
Moreover, the fixed length orphan spins n2 couple lin- 
early to the pp and thus appear as sources in this effec- 
tive action. Thus, we obtain 



Z cS cx JvnJvxJ Vpcxp + i ^ (J3h a + iA")(M~ 1 ) ff) ' (ph a + iAZ ) - i ^ p a ?g 



(11) 



Here the matrix M is defined by rewriting the quadratic 
part of iS e ff as 

YtftM^/l, (12) 

rr 1 a 

(3 is the inverse temperature, 

A2 = <S? A A£ + 5?f B t# o , (13) 

and h a may be assumed to be of the form h a = h5 az 
without loss of generality. 

Crucially, this may be rewritten in terms of the suscep- 
tibility matrix x whose elements Xff 1 (^) are the effective 
field theory prediction for the linear susceptibility of the 
spin at site r of the pure system in response to a local 



field applied at site r . To see this, we first note that 
A/ -1 can be identified with 

{ftf^^C&Sap, (14) 

the zero field correlations of the 4> field in the pure prob- 
lem: 

= C# (15) 
Furthermore, we have quite generally 

X* = PC* (16) 
and therefore, we may write 

(M"!) rV = Tx& (17) 
In other words, our effective action can be rewritten as 
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Z cff ex Jvnjvxj P^cxp j +^ J2 ( h ° + lTA ?)Xrr^ (h a + iTAZ) - i £ p%n% j 

\ rr a ? a J 

= jvnjv\f P/lcxp l^^x-^ A2C rV A^ +z^A2(^(^ x ^)-^ o <jj (18) 



where we have written the quadratic part of the action 
for A in terms of the correlation matrix C, while writing 
other terms using the susceptibility matrix x, since this 
proves to be the most physically transparent representa- 
tion for the subsequent analysis. 

One may now integrate over the Lagrange multipliers 
A2 and /i2 to obtain an effective action that couples all 
the orphan spins Sn~ to each other and to the external 
magnetic field h. The structure of the resulting effec- 
tive action for the n is controlled by the structure of the 
quadratic form T'\CV\ 1 where the projection operator 
Vk restricts C to the subspace spanned by r for which 
A^ is non-zero, i.e r corresponding to orphan sites r a and 
vacancy sites r v . 

From knowledge of this quadratic form, which depends 
crucially on the disorder configuration, one can compute 
(n~ ){h,T) of each orphan spin at a given temperature 
and field. Finally, by analyzing the path integral over 
</> with the orphan spins replaced by sources of strength 
equal to the orphan spin polarization (np ){h,T), we ob- 
tain {<pp){h 1 T), the spin texture induced around itself by 
each orphan. 

In Section I VII we summarize the results of such an 
analysis for one, two and three orphan configurations, 

I 

Z c ff oc J dnS(n 2 — 1) J dflcxp ^~~2 



relegating details to Section IIVI and Appendix [5] 

IV. CALCULATIONS WITHIN EFFECTIVE 
FIELD THEORY 

Since the effective field theory approach outlined above 
is somewhat novel, it is perhaps useful to first address 
some "structural" questions regarding the overall logic 
of this approach and its consistency before we turn to 
the actual calculations. 



A. Treating one spin as a unit vector in absence of 
orphans 

To this end, we first consider the pure system and ask: 
What would this approach predict if we chose to retain 
any one spin, say the one at site r, as a fixed-length vector 
n, while using the effective field theory description for the 
rest of the system? To answer this, one simply notes that 
Cp r = 1/3 independent of T, and that j xpr = Xf(T), 
the susceptibility per site obtained within the effective 
field theory. Xr{T) equals Xkag(T) (susceptibility of a 
site in the kagome layers) or Xap(T) (susceptibility of a 
site in the apical layer) depending on the location of r. 
Our approach would give 



rr a / 

(19) 



Note that the 0(n 2 ) terms add up to a constant because 
of the unit length constraint on n. From the above equa- 
tion, we get that (n z ) = (^-^ze + coth(3Xf?/i)^ which re- 
duces to (n z ) = x?h in the low-field limit. 

Another natural question that comes to mind is the 
following: What does this approach predict for an arbi- 
trary site r if we choose to retain the spin at that site as 
a fixed length vector n in a sample with a single vacancy 
at the origin? VaCVa is now two dimensional (corre- 
sponding to the vacancy site and the site r), and we 



obtain the following effective action for the spin at r: 

Z off oc J dnS(n 2 - 1) exp(C c ff (r, 0)h ■ ft) (20) 

where C cff (r,0) = ( X r - 3C fXo)/(§ - 3C&,). Crucially 
C e ff (r, 0) stays finite for all r as T — > ruling out Curie 
response of any kind in the single vacancy case. When r 
is well seperated from 0, we can ignore Cqp and C e ff = 
3xpi which is identical to the result for the undiluted 
problem derived above. When both sites r and belong 
to one of the kagome layers, C c ff gets further simplicd 
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to C e ff = i+zc,^ - Thus the local magnetization (n z (r)) 
at low T and uniform field h encodes the information 
about the spin-spin correlations of the parent spin liquid 
(i.e. Co?) in the case of a single vacancy. However, the 
orphan spins lead to parametrically stronger effects and 
the effects of single vacancies in a diluted lattice can be 
ignored to a good approximation. 

In both these examples (no vacancy and a single va- 
cancy), treating one spin as a fixed-length vector while 
describing the rest of the system by a Gaussian field re- 
sulted in a prediction that the spin that was singled out 
continues to have a finite magnetic suscepbility in the low 
temperature limit. When would our approach predict a 
Curie-like response with a divergent susceptibility in the 
low temperature limit? 

The answer has to do with the low temperature limit 
of the eigenspectrum of VaCVa- When the projector Va 
projects onto the subspace spanned by all sites of a sim- 
plex (i.e three sites that form a triangular simplex, or 
four sites that form a tetrahedral simplex), it is easy to 
see that the resulting VaQVk has one eigenvalue which 
goes to zero with temperature. The corresponding eigen- 
vector is the uniform eigenvector with equal amplitude 
at all sites of the simplex — this is in effect a consequence 

of the fact that the Hamiltonian % ([4>rfj enforces the 

constraints 

£ #=0 and ^# = (21) 

which in turn guarantees the fact that the vector with 
equal amplitude at all sites of a simplex is a zero mode 

I 



where both C' nn (nearest neighbor correlation on the 
kagome layer) and ((0^) 2 ) are calculated in the parent 
spin liquid state at h = 0. At low T, the above effective 



which immediately shows that removing two (but not 
one) spins from the triangular simplex leads to the ab- 
sence of all the pair-wise interaction terms above and the 
surviving h ■ n term then leads to the Curie response of 
the orphan spin. Also note that at low temperatures, 
the orphan spin acts like a free spin in a magnetic field 
of "ft,/2" instead of the applied field h. This is a direct 



of VaGVa- At finite temperature T, the precise state- 
ment is that the eigenvalue eo corresponding to this mode 
scales with temperature as eo = T. 

When there is more than one orphan spin present in 
the system, every defective simplex gives rise to a Curie- 
tail in the spin susceptibility even in the presence of the 
other defective simplices — the correlations between dif- 
ferent defective simplices will show up as entropically 
generated effective interactions between the orphan spins 
on these defective simplices; these are the effective ex- 
change couplings J e ff that we calculate later in this ar- 
ticle. Obtaining the low temperature behaviour of the 
orphan spins SH? g from knowledge of these interactions, 
one can in principle feed this information back in to com- 
pute the expected response of the surrounding spin liq- 
uid, and thus obtain the low-temperature properties of 
this system of interacting orphan-texture complexes. 



B. A single-simplex embedded in the spin-liquid 

Another way to re-phrase the argument is to consider 
three unit-length vectors on the three sites of a triangular 
simplex. We can then integrate out the rest of the soft 
spin degrees of freedom to obtain an effective action for 
these three spins on the simplex. This is then the "single- 
unit" action in this effective field theory description: 



(22) 

I 

action can be simplied by noting that ((</>^) 2 ) = T by 
equipartition, C nn = -± + 0(T) and Xkag = | + 0(T). 
Then, we get the following action 



(23) 

I 

consequence of the rest of the field being screened by the 
coupling to the surrounding spin liquid. 



Z e g oc / drTidnvdria exp 



mF)h . {nl+n2+n , )+ -—^ 



Hi ■ n 2 + n 2 ■ n 3 + n 3 • Hi 



Z cS oc / dnidn 2 dn 3 exp ( -(3h ■ (Hi + n 2 + n 3 ) - (3(ni ■ n 2 + n 2 ■ n 3 + ri 3 ■ Hi) 
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C. Texture induced by orphan spin 

The induced spin texture at a point r 2 far away from 
an orphan spin located at r\ can be expressed simply 
in terms of the parent spin liquid properties: When 
V2 — r± I S> a, where a is the lattice spacing, the detailed 
"internal" structure of the orphan spin simplex becomes 



irrelevant and we can simply impose </>a(^i) = ^orphan 
instead of removing two vacancies from the simplex. To 
calculate the spin polarization at r 2 , it is useful to first 
integrate out all other (r ^ r^) <fi(r) to obtain an effec- 
tive action that couples </>(r 2 ) to the unit vector n orp h a n- 
In practice, we do this by integrating over all 4>{r) while 
constraining </>(r 2 ) to take on a fixed value. This gives 



Vn 



orphan 



r,r — ri,r2 



(24) 



where /i2 at r ~ r\ , r 2 respectively impose the constraints 
that (f>A on the "orphan" triangle at fi equals n orp h a n and 
0(r 2 ) is held fixed. In the above, the matrix 

/ (<M(ri)-0(f 2 )) \ 

\ <0A(fl)-0(f2)) 1 / 

V 3 3/1 



is our approximation to P\CPa obtained by ignoring the 
internal structure of the orphan spin simplex (correlators 
appearing in the matrix elements represent correlations 
of the pure spin effective field theory in zero field). Per- 
forming the integrals over the /2 fields, we finally obtain 



ry f ~>-. r f ^\ / 3(</>A(fl) • </>(f 2 )) - Xl^\ I ^ ka 9t - (~>\ 

Z ef f oc / X>ri orphan (ri)I?0(r2) cxp — - — — ■n orphan (r 1 ) • </>(r 2 ) + ^ h ■ n^^irx) (25) 

J V^A> - (0A(ri) -0(r 2 )) 2 (<^) 

2((^)-(^(r 1 )-^ 2 ))^) W2j 



We now compute both (^o rp h an (^i)) an0 - {^{^2)) from 
this effective action. At low temperature, we obtain 
( n orph an (^i)) = B{h/2,T), the polarization of a unit- 
length spin at temperature T in response to an external 
field of magnitude h/2; this was also obtained from the 
more detailed calculation in the previous section where 
the two vacancies on a triangular simplex were explictly 
incorporated to compute the orphan spin response. At 
vanishingly small fields h deep in the low temperature 
spin-liquid regime (h <C T <C 1), we obtain: 



(4> z {r2)) «x(r 2 )/i + 



T 



B(h/2,T) .(26) 



Thus, we see that the texture induced around an or- 
phan spin is intimately related to the spin- "charge" cor- 
relation function of the parent spin liquid (where 0a is 
the thermally generated vector- "charge" defined earlier) . 



In the Appendix, we use a general long-wavelength anal- 
ysis of the properties of a "Coulomb-liquid" with fluctu- 
ating fields to gain insight into the nature of these spin- 
"charge" correlations. 

D. Orphan spin interaction 



As the orphan spin textures are extended degrees of 
freedom, it is a priori not at all obvious how they will 
interact as the centers of two of them are brought closer 
together. The great advantage of the method we have 
developed above is that it readily generalizes to the case 
of any number of orphan spins, albeit with increased cal- 
culational effort. We document the details relevant to 
the two orphan case in Appendix. [C] The final form of 
the action is about as simple as could have been hoped 
for: 



/ T>n\Dn2 cxp (- 



■(3J e ffHi ■ n,2 + fihs\n\ + /3/is 2 n 2 ) 



(27) 



Basically, the two Zeeman terms for the individual or- mented by an effective exchange 'constant' J e ff, which 
phan spins, where Si,s 2 — ¥ 1/2 at low T are supple- 
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depends on the location of the orphans, and on temper- 
ature. 

l3J e ff{f, T) can be appoximately expressed in an illu- 
minating form, which becomes exact in the limit r a, 
i.e. when the two orphan sites are well seperated. When 
r > a, we ignore the "internal" structure of the or- 
phans as before and simply impose <Pa(^i) = n\ an d 



</>a(^2) = ^2 instead of removing two spins from each 
of the triangular simplices on which the orphans reside. 
Then integrating out the remaining degrees of freedom 
yields the following effective action (here we stick to 
h = for notational simplicity): 



J VnVpS(nl- 1)6(4 -l) ex p l~ E E ^" ( 28 ) 



r,r — ri,r2 



where the Lagrange multipliers p impose the constraint 
that the total vector spin on the two "orphan" triangles 
equals Hi and n-i- The matrix M 2 above is a 2 x 2 matrix 
of the following form: 



appropriate linear combinations of the spin-spin correla- 
tions. However, as we will see in the next section, their 
behaviour is quite different from the 1/r 2 behaviour of 
the spin correlations of the pure system. 



M 2 = 



(<fc(ri)) 

?a(t ? i)-0a(t ? 2)) 



where the correlators are again calculated in the parent 
spin liquid state in the absence of disorder (hence, the 
diagonal terms arc independent of r). Now, it is easy to 
integrate out the p, fields to obtain the effective interac- 
tion between the orphans: 



J eff 



OC 



T>n{Dn2 exp (— /3 J e jjfi\ ■ n%) (29) 



where 



eff 



i>A(ri) ■ <M(r 2 )) 



(30) 



Thus, PJ e ff is determined by the "charge" - "charge" cor- 
relator of the parent spin liquid. In the Appendix, we use 
a general long-wavelength analysis of the properties of a 
"Coulomb-liquid" with fluctuating fields to gain insight 
into the nature of these "charge" - "charge" correlations. 

To summarize this section, the effects of putting non- 
magnetic impurities in the parent spin liquid show up in 
the following manner: Single vacancies leads to a spin 
texture that follows the spin-spin correlation (Eq |20|) of 
the parent spin liquid. Orphan spins, generated when 
all but one spin are substituted by non-magnetic impu- 
rities in a simplex, however generate a texture that de- 
pends on the spin- "charge" correlation (Eq [26)) of the 
parent spin liquid, where the "charge" is located on the 
orphan spin simplex. The orphan spins interact via a 
pair- wise Heisenberg interaction J e f f that is essentially 
determined by the "charge" - "charge" correlator (Eq [30]) 
of the parent spin liquid, in which the two "charges" are 
located on the two orphan simplices. Note that the spin- 
"chargc" and "charge" -"charge" correlations are simply 



RESULTS ON THE PURE SYSTEM 



Thermodynamics: Field theory and 
Monte-Carlo simulations 



The effective field theory results can be worked out 
for a pure system as the lattice is sufficiently symmetric 
to permit an analytical treatment despite its non-trivial 
seven-site basis. Working in Fourier space with this seven 
site basis, we obtain the following expressions for the 
magnetization of the kagome layer spins mkag, and the 
apical spins m ap : 



m kag = - 



m 



ap 




(31) 



Since p\ , P2 > at any T, this immediately implies that 
mkag,m ap > for all T within the largc-N approxima- 
tion. At T = 0, mkag = § and m ap = 0. The leading 
temperature corrections to this T = result are 



mka 



T 
2" 



P2 



lap = ~ (PiT H ) 



(32) 



Upon comparing with MC results for 0(3) spins (see 
Fig.|4]and Fig. [5]), we find that these effective field theory 
results for mkag and m ap exhibit important qualitative 
departure from the actual behaviour of 0(3) spins. First, 
from a comparison with the MC data, we see that the 
field theory prediction for mkag has the right dependence 
on T but with the wrong sign of the coefficient of the lead- 
ing low-temperature correction. Second wc note that our 



MC results demonstrate that low temperature m ap rep- 
resents a diamagnetic response to an external field, in 
the sense that the magnetization develops in a direction 
antiparallel to the applied field. This diamagnetic correc- 
tion is seen to be non-analytic, m ap ~ —\T, and cannot 
be captured within our effective field theory. 

To understand these failures of the effective field theory 
better, we have also studied a classical 0(4) Hcisenberg 
model on the same lattice using MC methods. From our 
MC results for the 0(4) case, we see that Xkag indeed 
follows the effective field theory temperature correction 
in Eq[32]at low temperature. Also, the diamagmetism of 
the apical layer spins is much reduced and follows m ap ~ 
— T at very low T. 

We would like to stress that a subleading correction 
due to thermal order-by-disorder effects is in keeping with 
the expectation from constraint-counting which is for the 
paramagnetic regime to persist all the way to T = in 
SCG0<2i In particular, one expects SCGO 0(3) spins 
to exhibit a spin liquid phase with a low-temperature 
specific heat of 6/7 per spin, indicative of two zero-energy 
modes per unit cell. This prediction is in quantitative 
agreement with our Monte Carlo results, Fig. [51 which in 
addition shows no sign of a phase transition down to the 
very lowest accessed temperature. 



B. Dynamics: Monte-Carlo simulations 

A conceptually connected but nevertheless distinct di- 
agnostic for the spin liquid regime at low temperature 
consists of considering the dynamics of the syste m 21 ' 22 , 
which greatly differs between frustrated and unfrustratcd 
systems 2 ^. The simple prediction is for the autocorrela- 
tion function of the spins to decay exponentially, on a 
timcscale set by the inverse temperature. Whereas we 
have not done molecular dynamics simulations of the 
proper equations of motion including all the relevant con- 
servation laws, our Monte Carlo results, obtained in sim- 
ulations that only use strictly local single-spin-flip moves 
(i.e., with the "microcanonical" over-relaxation moves 
switched off), are nonetheless strong evidence that the 
system does not enter a glassy state. Instead, the au- 
tocorrelation function Ck ag (t) = ^J2p(Sr(t) ■ SV(O)), 
where the sum is over spins on the two Kagome layers 
exhibits an exponential decay at a timescale which grows 
algebraically with the inverse temperature (Fig. [3]), in 
fact being precisely proportional to the inverse temper- 
ature as expected in a spin-liquid phase. Similar results 
were obtained for the apical spins (not shown). 

All of this taken together implies that SCGO is a model 
system with a spin liquid phase closely related to that of 
classical Hcisenberg magnets on the pyrochlore lattice, 
and other magnets in which thermal order-by-disorder 
effects do not produce an ordered state. Furthermore, as 
we will see below, the effects of vacancies are in several 
ways stronger in d = 2 than in d = 3; this, together with 
detailed experimental results available in the literature, 




T 



FIG. 2. Specific heat C v shows no signs of a phase transition 
for x — pure SCGO lattice as the temperature T is lowered. 
It converges to the value 6/7 obtained from a mode-counting 



1.0 


1 1 1 




\ — (5=4096 




- \. — (3=2048 


0.8 


— (3=1024 




— (3=512 






^0.6 








u 




0.4 




0.2 


1,1 



I i I i L 

2 4 

Tt 



FIG. 3. Autocorrelation function Cfc a9 (t) = (S(t)-S(0)) where 
S is a spin in the kagome layers for a system of size L = 32 at 
x = 0. Note that the autocorrelation time scales as ^ at low 
T, which shows that the x — system remains a spin liquid 
down to the lowest temperatures accessed. 

provides much of the motivation of studying SCGO as a 
candidate spin liquid with quenched disorder. 



VI. RESULTS ON IMPURITY EFFECTS 

A. A pair of vacancies on the same triangle: Single 
orphan physics 

To obtain the behavior of a single orphan texture, we 
start with the pure system and remove two magnetic sites 
and use the general procedure discussed earlier to obtain 
the effective action of the orphan spin fEq |2"2"j) . From this, 
we obtain that the orphan spin acts like a free spin a field 
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FIG. 4. Zero field susceptibility of the kagome layer spins (top 
panel) and the apical spins (bottom panel) for 0(3) and 0(4) 
spins. The apical spins show a —aVT behavior at low T for 
the case of 0(3) spins. 
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FIG. 5. Total susceptibility for 0(3) and 0(4) spins. Note 
that the results for 0(4) spins are in qualitative agreement 
with expectations from large- TV theory. 



h/2 at low T and hence, (n z orphan ) = B(h/2, T). We then 
use this result to calculate the full texture on the lattice 
scale using the approach detailed in the previous section 
and appendix. The screening of half of the magnetic field 
for the orphan spin at low T and the detailed form of 
the texture obtained from this procedure both agree very 
well with the MC results for the 0(3) model with two 
vacancies on a triangular simplex; a summary of these 
results has already appeared in our Letter—. 

Here we explore the result further by connecting it to 
correlations of the pure system. We cast the solution for 
the orphan induced spin texture in a fairly compact form 
using appropriate correlation functions of the parent spin 
liquid: 



<(0a) 2 > 



6(1/3 - C nn ) 
3CL 



1/3 - C n 



B(h/2,T) - 3h Xkag 

(33) 



where (</3 z {f\)4> z {?%)) is the spin-spin correlation between 
the spin at r\ (position of the orphan spin) and the spin 
at r*2 in the parent spin liquid and (<f> z ^(fi)(f) z (r2)) is the 
spin- "charge" correlation where <j>%{r\) is the z compo- 
nent of the sum of the three soft-spins on the orphan sim- 
plex (in the parent spin liquid) . At low T and sufficiently 
far away from the orphan spin, the above expression can 
be further simplified to give 

(^(ft))*x(5)ft+ — y " B{h/2,T) (34) 

This latter form can also be derived directly by argu- 
ing that the internal structure of the defective simplex is 
unimportant for r ^S> a (Sec IIVI) . 

The comparison between the approximate answer (Eq 
|3"4"1) and the full effective field theory result (Eq |3"3"|) is 
shown for a system of size L — 50 in Fig [6j As is clear 
from the Figure, the approximation only deviates signifi- 
cantly from the full result when the texture very close to 
the orphan spin is considered. Further, for |fi — fa| S> a 
(where a is the lattice spacing), the spin- "charge" cor- 
relator (</>2,(fi)0"(r2)) is expected to satisfy the scaling 
form 

(^(fi)^(f 2 )) = vin^F.dn - r 2 )VT) , (35) 

where F\(x) ~ l/l^l when |a?| <C 1 and decays exponen- 
tially at large \x\, and r\ is the sublattice index of the 
bipartite dual lattice of simplices, taking on a value +1 
for the A-sublattice, and —1 for the B-sublattice; the 
rationale behind this expectation is detailed in the ap- 
pendix and relies on our analysis of the long-wavelength 
properties of a "Coulomb liquid" with fluctuating fields. 
In Fig we see that this expectation is borne out by our 
results. 
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FIG. 6. The texture obtained from the spin- "charge" corre- 
lator of the parent spin liquid in the Kagome bilayer shown 
along a one-dimensional cut on the lattice near the orphan 
spin (see inset). The "charge" location is fixed on a triangu- 
lar simplex as shown in the inset. The approximate version 
of the effective field theory result only deviates significantly 
from the full prediction (labeled "large- iV" in figure) close to 
the orphan spin site. System size is L = 50 at f3J — 2048 and 
h/J = 0.125. 



B. Form of the screening cloud 




FIG. 7. The spin- "charge" correlator within our effective field 
theory is proportional to T 3 ^ 2 Fi(f\/T), as expected from the 
Coulomb liquid analogy (see Appendix). 

over which the total spin of the zero temperature orphan- 
texture complex reaches S/2 to the same accuracy. The 
robustness seen at finite temperature is thus connected 
with the moment of the zero temperature orphan-texture 
complex approaching its asymptotic value rather quickly 
in the sense of Eqn. as one goes further and further 
out from the core of this complex. 



An interesting aspect of our results is the fact that the 
texture and the resulting impurity susceptibility deviates 
so little from the asymptotic low temperature predictions 
even at sizeable temperatures of order 0.1JS 2 ; this was 
already noted in our earlier Letter—. To understand this 
robustness better, it is useful to consider just how the 
orphan-texture complex acquires a net spin of Sf ot = S/2 
at T = as more and more spins around the orphan are 
taken into account. 

To explore this, we begin by noting that the orphan 
spin is fully polarized by an infinitesimal magnetic field 
at T = 0, and produces a staggered spin texture around 
it that decays as 1/r from the Coulomb phase anal- 
ogy. Defining a smeared total spin operator S' ot (^) = 
^2 i rS z (r) exp(— r 2 /£ 2 ), where r is measured from the or- 
phan spin site, we find that 

sU0 = ^ + Sf(0 

where /(£) ~ V£ 2 , £ > 1 (36) 

for a two-dimensional Coulomb spin liquid (sec Appendix 
for the calculation of this operator in the simpler case of 
the planar pyrochlore lattice). Thus, S* ot (£) approaches 
S/2 quite quickly (£ ~ 6 is enough to approach within 
1% for the planar pyrochlore lattice). 

This has important implications for the finite tem- 
perature properties. The magnetic susceptibility can be 
though of as arising from spin-5'/2 orphan-texture com- 
plexes to a given accuracy as long as the thermal correla- 
tion length £(T) ~ \j\fT are larger than the length scale 



C. Two orphans: Effective interactions between 
orphan spin textures 

The interactions between the orphan spins can be cal- 
culated by considering the two orphans as fixed-length 
vectors and integrating out the rest of the soft-spin de- 
grees of freedom <p(f) as detailed in the previous section 
and appendix. In this way, we obtain 

J eff « ZI^MJ^M (37 ) 

Further, for |r*i — r^l 3> a, where a is the lattice spac- 
ing, we expect that the "charge" - "charge" correlator 
(</>2(ri)<^(r 2 )) satisfies a scaling form 

(^(ri)^(f* 2 )) = -rtnW^Fdn-r^VTps) 

at low temperatures, where rj is the sublattice index of 
the bipartite dual lattice of simplices, taking on a value 
+1 for the A-sublatticc, and —1 for the £?-sublattice; 
the rationale behind this expectation is again detailed 
in the appendix and relies on our analysis of the long- 
wavelength properties of a "Coulomb liquid" with fluc- 
tuating fields. From Fig [31 we see that this expectation 
is borne out by our results. Since {<p\) = 3T at low T, 
this implies a scaling form 

Jeffin - r 2 ,T) = rtnW^TJiVTin - r 2 )) (39) 
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FIG. 8. (Top panel) "Charge" correlator calculated between 
triangular simplices on the same layer in the SCGO lattice. 
(Bottom panel) The effective interaction for two orphans 
in the same layer (upper curve) and different layers (lower 
curve). System size used is L = 200. 



textures themselves, which are screened at T = when 
the orphans are in opposite layers, since a "charge" zero 
"dipole" formed by two orphans on opposite layers leads 
to a 1/r 2 far-field behavior instead of the 1 jr profile of a 
single spin texture at T = 0. 

Fig. [5] shows that the effective field theory computa- 
tion for J e f f (r , T) agrees very well with the effective in- 
teraction obtained from direct simulations of the 0(3) 
Heisenberg model with vacancies. In the simulations, 
we create two orphan spins by removing two vacancies 
each from the chosen triangular simplices. Because of 
the complicated geometry of the lattice, many symmetry 
inequivalent choices are possible and here we show three 
of them. We monitor (S(Q) ■ S(r)) in the simulations at 
different r and various (low) temperatures, where S(0) 
and S(r) refer to the two orphan spins. This quantity is 
then computed using the effective field theory result for 
J e ff(f, T) and the agreement is excellent in all the cases. 
The effective field theory computations which were done 
on finite lattices for SCGO, also capture the finite-size 
effects in the system very well. 

The probability distribution of a; = S(Q)-S(r) obtained 
from the MC simulations (see Fig ITTJ|) can also be fully 
matched to P(x) oc exp(— J e fj(f, T)x) to rule out inter- 
action terms of the form (Si • S2) 2 which are not for- 
bidden on symmetry grounds, but appear to be absent. 
Finally, we note that one may in principle plug this infor- 
mation back in and obtain the response of the surround- 
ing spin liquid to this pair of interacting orphan spins, 
and thereby compute the low-temperature behaviour of 
this system of two interacting orphan-texture complexes 
(as emphasized earlier in our detailed summary of the 
effective field theory computations). 



where J{y) is seen to have the asymptotic behaviour 

J{y) ~ log \y\ for \y\ < 1 

J(y) ~ exp(-|ifl) for \y\ > 1 (40) 

The dependence on sublattice index 77 leads to another 
interesting observation upon noting that all triangular 
simplices on the upper Kagome layer have ?; = +1, while 
all triangular simplices on the lower Kagome layer have 
i] = —1. Therefore, two orphans in the same layer in- 
teract antifcrromagnctically and have a vanishing net 
Curie response in the limit of low fields and temperatures 
smaller than this interaction scale. On the other hand, 
two orphans in opposite layers couple ferromagnetically, 
leading to an enhanced Curie tail due to a 'restituted' 
moment equal to that of a full free spin SI Note that 
this restituted moment arises because most of the spin- 
density that leads to the fractional moment of S/2 for a 
single orphan is localized close to it, as we discussed in 
the previous section. It is also interesting to note that 
this behaviour is in sharp contrast to that of the spin 



D. Three orphans: Absence of multi-spin 
interactions 

From the structure of the effective action for the 
Lagrange-multiplier fields detailed earlier, it is clear that 
our effective field theory always gives pair-wise interac- 
tions even when the number of orphan spins is greater 
than two. This is a strong prediction. The way we 
check this from our numerics is to place the three orphan 
spins in a symmetric equilateral triangle configuration as 
shown in Fig 111! 

The first calculation we do is to calculate ((Si + S2 + 
S3) 2 ) for a given system size L and inverse temperature 
(3 from the effective field theory and check it with the 
result obtained from a pair- wise J e //(f) interaction (see 
previous section). The agreement is extremely good. We 
show in the table below the results from runs at L = 32 
at f3J — 256 for four different separations (see Fig II II for 
the configuration chosen in these runs). 
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FIG. 10. Probability distribution of the dot product of two 
orphan spins measured in MC simulations (points) exactly 
matches the result obtained (solid lines) by using a simple 
Heisenberg interaction term having magnitude and sign pre- 
dicted by the effective field theory. 




FIG. 11. Three orphan spins placed in an equilateral triangle 
arrangement. The circles indicate the vacancies. 



FIG. 9. Agreement between effective field theory predictions 
(solid lines) for orphan spin correlators, and actual results 
(points with error bars) for the same quantities obtained from 
MC for the 0(3) system shown for three inequivalent orphan 
spin placements (shown in corresponding insets). 
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Secondly, we probe the probability distribution of the 
orphan spins to see how well a pair- wise interaction pic- 
ture can explain it. Since, the full probability distribu- 
tion function is complicated even for pair-wise interac- 
tions, we fix two of the orphan spins 5*2 — S3 = z in 
the MC numerics and then monitor the probability dis- 
tribution Sf of the third orphan spin. If the pair-wise 
interaction picture is true, then the resulting distribution 
will be proportional to exp(— 2J e f f(r)S[), where J e //(^) 
is the pair-wise interaction strength, and this is exactly 
what we observe from the simulations (Fig IT2"1) . 
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FIG. 12. Probability distribution of Sf in the three orphan 
problem when the other two orphans are fixed in the z direc- 
tion (as detailed in text). Points represent values measured in 
MC simulations, while lines are fits obtained from the effective 
field theory prediction of purely bilinear two-body exchange 
interactions. 

VII. GENERALISATIONS 
A. Orphan tetrahedra 

The above arguments all hold for orphan (thrice- 
defective) tetrahedra as well as the orphan (twice- 
defective) triangles discussed throughout. The emergent 
gauge charge of an orphan tetrahedron is opposite to that 
of a triangle in the same layer. In a random dilution 
model, the probability of the former is 3a; 2 (1 — x) while 
that of the latter is 4x 3 (l — x), much smaller in the limit 
of small x. 



B. Other lattices, and dimensions 

Several central results readily generalise to other lat- 
tices: orphan spins can be induced in 0(n > 4) spin 
liquids on the kagome lattice, and for 0(n > 3) spin liq- 
uids on the pyrochlore lattice, which is of course where 
they were first identifie d 15 ' 16 . Carrying emergent gauge 
charges, they interact via a r~ d+2 Coulomb interaction 
for d > 3, the generalisation of the d = 2 logarithm 

implied by Eq. 1401 analogously, their textures decay as 
r -d+i_ 



VIII. CONCLUSIONS AND OUTLOOK 

Motivated by the extensive set of experimental data 
available for SCGO, we have studied the Heisenberg 
model on the corresponding lattice as a model system for 
a spin liquid. We have established in detail that SCGO 
remains in a spin liquid phase down to the lowest tem- 



peratures, as expectation based from mode counting ar- 
guments. 

We have then studied in detail the response of this spin 
liquid to the inclusion of disorder in the form of vacancies. 
By developing a field theory which captures the hard-spin 
nature of spins near a vacancy and treats entropic effects 
on the same footing as energetics, we were able to get 
an analytical handle on the resulting phenomena — these 
predictions were found to be in excellent agreement with 
direct Monte-Carlo simulations of the Heisenberg model 
with vacancies. 

In particular, we have found that vacancies that leave 
behind more than two spins in a simplex generically lead 
to a regular low-temperature limit of the susceptibility. 
On the other hand, the presence of an orphan spin, all of 
whose neighbouring spins in a simplex have been removed 
by dilution, leads to a Curie tail in the susceptibility, with 
a characteristic 1/T divergence in the low temperature 
limit, corresponding to the susceptibility of a free "spin 
S/2". This fractional moment occurs as a combination 
of two effects: First, the coupling to the surrounding 
spin liquid "screens out" half of the external field h seen 
by the orphan spin, so that it behaves as a spin S in 
a field h/2. And second, the surrounding spin texture 
develops a diamagnetic response that "cancels off" half 
the polarization of the orphan spin. 

These orphan-texture complexes experience long-range 
mutual interactions on account of their extended nature, 
captured intuitively by an analogy to the electrostatics 
of the Coulomb phase of the field theory describing the 
spin liquid. One of the central advances in this work is 
our derivation of the relevant scaling functions fully de- 
scribing vacancies in SCGO — these quantitatively cap- 
ture both thermal and energetic effects on an equal foot- 
ing, requiring only knowledge of the correlators of the 
pure system! 

The field of vacancies in unconventional magnetic 
states is a rich one with a long and interesting history, 
dating back (at least) all the way back to Villain's work 
on canted spin states in his seminal paper on insulating 
spin glasses 2 ^. This remains an exciting frontier, with 
many interesting avenues worth exploring. The obvious 
next step would focus on the many-body states result- 
ing from the orphan interactions described here. Under- 
standing this many-body physics is of considerable gen- 
eral interest, since the form of our interactions is actu- 
ally a quite general aspect of defects which cause viola- 
tions of the emergent Gauss law constraint in Coulomb 
phases. On general grounds, with such defects randomly 
distributed, this interaction can naturally lead to the ap- 
pearance of a spin-glass phase. We are currently inves- 
tigating this issue in detail^ 5 .. Returning to the specific 
case of SCGO, it will then remain to be seen whether 
the spin glass transition observed in experiments can be 
related to orphan-texture freezing. 

More broadly for the case of SCGO, we hope that our 
work will stimulate a more detailed study with the aim 
of better characterising the disorder present there, given 
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that our analysis of the NMR lincshapesi^ has found that 
the Curie tail cannot be explained with reference to or- 
phans induced by the nominal amount of uncorrclatcd 
vacancies in this series of compounds. In this context, it 
would be useful to characterize in more detail the statis- 
tics of substitution of Ga on the Cr sites, including possi- 
ble correlations, as well as better characterize other forms 
of disorder that may be playing an equally important role 
in the temperature regime studied experimentally. 
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Appendix A: Long- wavelength description of "charge" -"charge" and spin- "charge" correlators 

Our results rely crucially on the "charge" - "charge" and spin- "charge" correlators satisfying the scaling forms 

(^(n)0f(f 2 )) = rtir^T^F.dn - f 2 )VT) , 

(^(fi)0X(f 2 )) = -vinMr^Fttn - r 2 )VT) (Al) 

where a denotes the spin space index of 0, and by attaching a real-space index i, we have emphasized that (f> is best 
thought of as a spatial vector with orientation given by that of the corresponding dual lattice link on which this field 
lives. In the above, r] is the sublattice index of the bipartite dual lattice of simplices, taking on a value +1 for the 
^4-sublattice, and —1 for the B-sublattice. 

These scaling forms are a consequence of the fact that <p a behaves at low temperature like a fluctuating magnetic 
field on the links of the dual lattice, while t/*^ is proportional to the divergence of this fluctuating magnetic field, 
i.e. the fluctuating magnetic "charge" on a simplex . To understand this scaling behaviour, it is useful to consider 
a coarse-grained theory formulated directly in the continuum. The success of this continuum approach relies on the 
fact that the geometric details of the lattice only affect the short-distance form of these correlators (at scales r ~ a, 
where a is the lattice spacing) . The continuum theory detailed below is therefore expected to apply in a broad regime 
\f\/a ^> 1 with no restrictions on ^/^(T 1 ) ~ H'v/T'. 

With this in mind, we start with the continuum effective action for a fluctuating classical magnetic field 

s = \ J A dd < ga ) 2 + ^rJ A ddx{y ■ z a)2 (A2) 

where the repeated spin-space vector index a is summed over, and the subscript A serves to remind us that this theory 
is defined with an upper-cutoff A ~ 1/a in momentum space. In the above, the first term is cntropic in origin (and 
hence, it tends to a finite value as T — > 0), while the second represents the Boltzmann weight for creating magnetic 
"charges", A c and k are phenomenological constants related to the energy barrier for producing "charges" and the 
"magnetic permeability" of the medium, and we identify 

€ ~ B? (A3) 

with sign chosen so that the magnetic field points from a A-sublatticc simplex to a B-sublatticc simplex in the 
microscopic version of the theory. 

In order to work with this continuum action, we decompose the magnetic field into a curl-free pure gradient part, 
and a gradient free purely rotational part 

B a ~ (z x V)a Q - Vu Q (A4) 

and perform the functional integrals over a and u. Since all correlations considered here are diagonal in the spin space 
index a, which plays no role below, we drop it in the rest of this discussion. 
We now have the correspondence 

<^(r) or 4>a{t) ~ ?](r)V ■ B ~ -r/(r)V 2 u (A5) 

where r) is +1 (—1) if 4>^(r) or <j>&(r) represents the total spin of a A-sublattice (B-sublattice) simplex at r; thus, 
lattice scale distinctions only enter our theory through rj{r). 
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For the "charge" - "charge" correlators, this immediately gives in d = 2 

(^(n)^(r 2 )) or (tKn^im ~ TnWnfa) f d2 ^^^ + \ T — ( A6 ) 

which, apart from an extremely short-ranged part [represented in our continuum approach as a contribution propor- 
tional to Sjj^r — r ) = J Td d qexp(iq ■ (r — r))], reduces to 

(<t>®(n)<t>®(? 2 )) or (^(fi)^(f 2 )) ~ -v(n)v(r2)T 2 T ^M^^-^VT) 

7 A c g 2 + re 

which clearly yields the scaling form mentioned above. Furthermore, from the structure of the integral in the above, 
it is clear that F(x) has precisely the small and large x behaviour described in the main text. 
For the "charge" -spin correlator, we have in d = 2 

J A c q^ + re 

which immediately implies the scaling behaviour quoted at the outset of this appendix. 



Appendix B: Details of the orphan spin texture calculations 



The matrix describing the action (see Eq 2]) for the pure system (including Lagrange multipliers and a shift in the 
zero of energy to set the ground state energy to zero) is Fourier transformed to read 



/ i + -£2- in 

1 T 2/3 J 2 V 1 

i(l + e -**») h 



M(k) 



\ 



Pi 

2/3J 






1 1 

2 2 



i(l + e 1 ^) 

l(l + e t(fc«-fev)) 



) 1 + 

/ x 2/3J 






1 

2 







1 + -£3- 1(1 

± ~ 2/3J 2 v x 

l(l + e ifc„) 1 

\{l + e lk *) 1(1 + e 4 ( fe -- fe «)) 



pi 

2,3,7 



l(l + e -«.) 
.(I _|_ g-'( fe «i- fe y)) 
1 



Pi 

2/3 J 



P2 / 

20J / 



See Fig [T] for the numbering convention chosen for the 7-point unit cell and x and y vectors for the Bravais lattice. 
The Fourier transform is defined as 4> a {r) = ^gexp(— ifc • r)(f) a (k) where a = 0, 1, • ■ ■ , 6 is the sublatticc index 

and these are given the coordinate r of the direct Bravais lattice point, and N d = L 2 is the number of sites in the 
Bravais lattice. We denote the unitary matrix that diagonalizcs M(k) by U(k): U(\ty M(k)E/(k) = D(k), where -D(k) 
is a diagonal matrix. Then, going to the variables cr a (k), where 4> a (k) = U a p(k)ap{k) diagonalizes the quadratic 
interaction matrix. 

For calculating the orphan spin texture, we impose three constraints in the soft-spin calculation as explained in 



Section UlI CI two for the two spins being removed, and one that enforces the condition that 



orphan! 



B{h/2,T) as 



a hard constraint. Furthurmore, since we are interested in calculating the full texture, we introduce a "source field" 
Ja{r) which couples linearly to <fi%(f). Going to momentum space, and changing to cr a (k) variables then leads to the 
following path integral (for the z component, parallel to the direction of the external uniform magnetic field h) 



Z = I DaK(k)Dai(k)Dv 3 D\±D\ 5 cxp{-f3jY,DUm<(k)) 2 + K(fc)) 2 ) + Ph^/W d £ Uf a {k = 0)<r* \k = 0) 

a,0 



k . a 



k.a 



(/J3^e,3(^o,fc,a) + A 4 f7e,4(n),fc,a) + X 5 U et5 (f ,k,a))a^(k) 



J2(j^(k)^(k) + Ji(ky a (k))^ 

^(fi 3 U 0t 3(fo,k,a) + \ Jj O A{r ,k,a) + \ 5 U , 5 (r a , k, a))a^(k) - i\fN~ d ^B{h/2,T)] 



(Bl) 



where we assumed that the orphan spin lives on sublattice 3. /J3 imposes the length constraint on the orphan 
spin and A4 and A5 impose <$> z a (r) = for the two spins removed from the orphan simplex. Furthermore, J a (k) = 
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U a p{k)Jp{k) and 

U et3 (f ,k,a) = U 3 R a (k) cos(k ■ r ) + U^k) sin(k ■ r ) 

U e ,4,(ro, k, a) = Uf a {k) cos(fc • r + k y ) + U( a (k) sin(fc • r + k y ) 

Ue,5(ro, k, a) = U^ a (k) cos(fc ■ r + k x ) + Ul a (k) sin(fc • r + k x ) 

UoM^o, k, a) = £/^(fc) sin(fc • r ) - Ul a (k) cos(fc • f ) 

U 0t i{?o, k, a) = U± a (k) sin(fc • r + k y ) - U[ a (k) cos(fc ■ f + k y ) 

Uo, 5 {ro, k, a) = U^ a (k) sin(fc • r + k x ) - U^ a (k) cos(fc • r + k x ) (B2) 
Solving the above path integral to obtain Z(J), the spin texture can then be calculated by evaluating 

(<r£(fc)) = Kmj. —- dZ{J ] 
J ^ a 2Z(J)dJ£(k) 

The field ((j> z a {r)) can then be obtained from (cr^(fc)) and (cr^(fc)) which leads to the result displayed in Eq 1551 



Appendix C: Details of field theoretic calculation of interactions between orphans 



Here, we show some essential steps needed to obtain the effective interaction between two orphan spins. As discussed 
in Section llVDl we keep the orphan spins as unit vectors and intergate the rest of the soft-spins. Let us consider the 
specific case when the orphan spins are both placed on the sublattice a = 3 (other cases can be similarly considered). 
Then we impose six constraints, four for the four spins being removed and two for pinning the two orphan spins to be 
n\ and n\. We will need to remember that the (nf ) 2 and (n,^) 2 terms in the action can be combined with similar x 
and y terms and are unimportant since these are unit vectors. The only terms that will be generated in the effective 
action will be of the form hsinf, hs 2 n\ and Jeffnln^ (the corresponding x and y terms will combine with this to 
give the J e //ni • 712)- The full path integral (for the z part) is of the following form: 



J Da^i^Daii^D^D^DXsD^DXsDXe cxp[-(3jJ2D aa (k)((a^(k)) 2 + (a^(fc)) 2 ) - - i\i±y[W d nl 

k,a 

ph^/WdYjUfuik = 0)a%(k = 0) +i^ 1 J2(UeMriXa)a^(k) + U 0tS (f 1 ,k,a)aUk)) 

a fi k.a 

i\ 2 ^(t> e ,4(ri, k, a)a*(k) + U oA (f u k, a)al(k)) + 1X3 ]T(J7 e ,5(n, k, a)a*(k) + U fi{fi, k, a)cr^(fc)) 



A.' . <:v 



k.a 



X^e,3(r 2 , k, a)<r*(k) + U 0i3 (f 2 ,k, a)ai(k)) + iX 5 ^{U eA {r 2 , k, a)<r*(k) + U oA (r 2 ,k, a)a T a (k)) 



k.a 



k.a 



iX 6 ^(t/ e> 5(r 2 , k, a)a*{k) + U , 5 (f 2 ,k, a)ai(k))} 



(CI) 



Integrating out the a 11 ' 1 fields from the above path integral gives the following: 



Z oc / DmDX 2 DX 3 Dfi 4 DX 5 DX 6 exp -N d X T MX + 



where the interaction matrix 2M is defined below: 
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FIG. 13. The planar pyrochlore lattice shown above. The lattice can be represented as a Bravais lattice with a two-point basis 
as shown. The lattice translation vectors x, y are also shown. Three spins have to be removed from a simplex (open circles) to 
produce an orphan spin. 



and the coefficients c a = hxkag (h-x(r) in general). The above matrix is the 6x6 matrix P\CP\ for the two orphan 
case. The correlators in the matrix M are calculated in the parent spin liquid state (i.e. no disorder) in a zero 
magnetic field. 

A numerically simple method to solve the remaining A and fi integrals is by going to the diagonal basis for M. 
Thus, we use (below we use the notation \i\ = Ai and /i4 = A4 for brevity) 

V T MV = M, VV T =1, X a = V a f,\p (C3) 

to obtain the following: 

Zcx J DX 1 DX 2 DX 3 DX 4 DX 5 DX 6 cxp y~N d J2 M aa X 2 a + i^/W d ^ f^c^ a j A Q j 

Z(XCXP 1 V -^ 4M m J (C4) 
From this, we finally get the form 

Z (x exp {—Jeffn 7 {ri2 + j3hs\n\ + f3hs2ri2) (C5) 
from which we can read off Jeff, s-y and S2 

eff ~ ^ 2M aa 

a 

1 \ - c a VipV a p 



Sl 



2^ % M 0P 



1 \ - c a V a pV 4 

S2= 2^^T^- (C6) 



Appendix D: T = calculation of smeared spin density of orphan spin texture 

The planar pyrochlore lattice, where the four spins in an elementary unit interact equally strongly with each other 
(Fig I13p. presents a particularly simple case in two dimensions where the effective field theory can be solved easily 
because the diagonalizing matrix U(k) is a 2 x 2 matrix (instead of being 7 x 7 as in SCGO). Therefore, we focus on 
this tractable case to illustrate the behaviour of the smeared spin density of a T = orphan spin texture. 

Smeared total spin operator: As we have noted in the main text, an orphan spin in a spin S magnet induces a 
scale- free 1 jr texture around it at T = in the presence of even an infinitesimal magnetic field. Formally adding up 
the spin polarization at all sites in the system then leads to S% ot = 5/2 via the argument outlined earlier. Strictly 
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FIG. 14. Exact calculation of the smeared total spin operator S^ t(£) for the planar pyrochlore lattice at T = 0; values on the 
y axis are measured in units of S. The solution shows that the spin density is mostly concentrated in the neighborhood of the 
orphan spin and decays rapidly away from it. 



speaking, S% ot is only conditionally convergent in two dimensions. However, our finite temperature results demonstrate 
quite clearly that assigning a total spin of S z ot = S/2 to the orphan spin texture is meaningful in that the impurity 
susceptibility due to a single defective simplex is asymptotically equal to the susceptibility of a free spin S/2 in the 
low temperature limit. Moreover, as was already highlighted in our previous Letter, this asymptotic low temperature 
result remains surprisingly accurate all the way up to temperatures of order T ~ 0.1J. 

It is therefore very interesting to see how the total spin converges to S/2 when spins only within a certain radius of 
the orphan spin (which drives the texture) are taken into account. To this end, we define the following smeared total 
spin operator 

SUO = X>*(f)) exp(-(f - f ) 2 /£ 2 ) (Dl) 

r 

where we use a Gaussian centered at ro, the orphan spin location, to regulate the conditionally convergent sum. Thus, 
the above smeared operator effectively ignores the contribution of the spins located at distances much greater than £ 
from the orphan spin — note that £ is not the thermal correlation length, but a free parameter here. It is clear that 
when ^ 0, S* ot (Z) -> S. 

By way of illustration, we calculate (S z (r)) and this smeared total spin (S^ ot (^)) for the planar pyrochlore lattice 
in the thermodynamic limit at T = using a lattice Green's function technique. At T = 0, we have XVeH = 
and the orphan spin S z (fo) = S when h — > 0. On the dual square lattice, finding the texture (S z (f)) is equivalent 
to determining the current on each bond of the lattice given that three bonds of the lattice are removed (the three 
vacancies around the orphan spin) and there is an input current of S at r (the polarized orphan spin). The condition 
J2pe^ SZ = translates to current conservation at each site on the dual lattice. The current profile can be calculated 
by evaluating the lattice Green's function of the square lattice with three removed links which determines the potential 
on each site. The current on each bond is then just the potential difference across the bond. To calculate the lattice 
Green's function, we use the fact that the Green's function of the square lattice with three links removed can be easily 
expressed in terms of the lattice Green's function of the full square lattice (see J. Cserti, D. Gyula and P. Attila, Am. 
J. Phys. 70 (2), 153 (2002)). 

In this way, we obtain the result: 

SUZ) = |+5/(0 (D2) 

where /(£) w 0.16/£ 2 when £ 3> 1; this behaviour is displayed in Fig[T3J Thus, >S t z ot (£) approaches S/2 quickly (Fig[T3|) 
and it is already within 1% of its £ — > oo value when £ ~ 6 for the planar pyrochlore lattice. In other words, the spin 
density, which peaks at the orphan spin, is concentrated close to it and falls rapidly as the distance from the orphan 
increases. And clearly, it is this rapid fall-off that is at the root of the surprising robustness of the orphan physics all 
the way up to temperatures of order T ~ O.IJS 2 . 
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